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Abstract 

We study properties of gravitational waves (GWs) from rotating core-collapse of a 15Mq star by 
performing three-dimensional general-relativistic hydrodynamic simulations with an approximate 
neutrino transport. By parametrically changing the precollapse angular momentum, we focus on 
the effects of rotation on the GW signatures in the early postbounce evolution. Regarding three- 
flavor neutrino transport, we solve the energy-averaged set of radiation energy and momentum 
based on the Thome's momentum formalism. In addition to the gravitational quadrupole radia- 
tion from matter motions, we take into account GWs from anisotropic neutrino emission. With 
these computations, our results present evidence that non-axisymmetric instabilities play an essen- 
tial role in determining the GW signatures in the rotating postbounce evolution. For our rapidly 
rotating models (albeit not rapid to trigger the classical bar- mode instability), we show that prec- 
ollapse density inhomogeneities give rise to millisecond variations in the waveforms, which is only 
clearly visible to an observer along the rotational axis. During prompt convection, we find that 
the waveforms show narrow-band and highly quasi-periodic signals which persist until the end of 
simulations. We point out that such feature reflects the growth of the one-armed spiral modes 
that develop under the influence of the standing-accretion-shock instability and the low-T/|M^| 
instability. The typical frequency of the quasi-periodic waveforms can be well explained by the 
propagating acoustic waves between the stalled shock and the rotating proto-neutron star surface. 
Although the GW signals exhibit strong variability between the two polarizations and different 
observer directions, they are within the realm of next generation detectors such as by KAGRA and 
Advanced LIGO, if the source with sufficient angular momentum is located in our Galaxy. 

PACS numbers: 04.25.D-, 04.30.-w, 95.85.Sz, 97.60.Bw 



I. INTRODUCTION 



The seminal paper [l| on the gravitational- wave emission from stellar core-collapse started 
with the following two sentences: "The current and near-future efforts (Weiss, 1979;Epstein, 
1979; Weber, 1979; Douglass and Braginsky, 1979)1^^] to detect gravitational waves have 
increased the demand for reliable theoretical calculations of gravitational radiation from all 
possible energetic astrophysical sources. Among all potential sources, supernovae, possessing 
highly asymmetric cores and being situated in our own galaxy, have been thought to be the 
most promising candidates, probably already within the realm of second generation gravita- 
tional wave detectors." It may be surprising that only by updatingthe references in the late 
1970's to the most recent ones (e.g., Advanced LIGO/VIRGO ytzj and KAGRA [sj for an 
observational side), above statements are true timelesslju, outlining the final goal toward 
which all the relevant studies regarding the gravitational-wave (GW) signatures of core- 
collapse supernovae (CCSNe) have been heading since then (see recent reviews in |9l. Il0|). 

Following [l|, most of the theoretical predictions have focused on the GW signals from 
rotational collapse of the supernova cores (e.g., 1 ll-l23| and references therein). In this 
context, rapid rotation of the precollapse core leads to significant rotational flattening of 
the collapsing and bouncing core, leading to a time-dependent quadrupole (or higher) GW 
emission. The waveforms of the bounce signals are categorized into three types, namely 
types I, II, and III. In general, a realistic nuclear equation of state (EOS) is stiff enough 
to prohibit type III waveforms. A combination of general relativity and electron capture 
suppresses multiple bounce in the type II waveforms. Therefore the type I waveforms are 
considered to be most generic, which exhibits a sharp negative spike at bounce followed by 



a subsequent ring-down phase [17|, |l8|, |20| 



After core-bounce, anisotropic matter motions due to convection 24j-|27l| and the Standing 



Accretion-Shock- Instability (SASI, 28l-l32| ) . and anisotropy in neutrino emission 27|, 



30, 



33, 



34| are expected to be primary GW sources with comparable amplitudes to the bounce sig- 



nals (see [9|, llO| for more details). In general, these postbounce GW signals should be more 
difficult to detect compared to the bounce signals, since they change much more stochas- 
tically with time as a result of chaotically growing convection and SASI in the non-linear 



hydrodynamics ( |27l . 



28 



30 



32 



35|). However an encouraging finding is recently reported 



^ although the most promising source for the first detection may be most likely for compact binary mergers.. 



by 3J] who made a detailed analysis on the GW signals from exploding two-dimensional 
(2D) models with currently best available neutrino transport scheme. They pointed out that 
an episode of convective activities in the vicinity of the proto-neutron star (PNS) |32| that 
trigger the onset of explosion is encoded in the GW spectrogram, which could be read out 
by a power-excess method |36l]. A combined analysis of the emergent neutrino signals [37 1 
should provide an important clue to the supernova engine in the context of the neutrino 
heating mechanism (e.g. 



38 



j|-l40| ) . which is otherwise very hard to access for conventional 
astronomy by electromagnetic waves (e.g., 41|). 

Regarding the bounce signalcl, rapid rotation is necessary for producing distinct GW 
bursts and it is likely to obtain ~ 1% of massive star population (e.g., 42|). Minor as they 
may be, a high angular momentum of the precollapse core is essential for the working of 
the magnetohydrodynamic (MHD) mechanism 43N48| . This is because the MHD mecha- 
nism relies on the extraction of rotational free energy of the collapsing core by means of the 
field- wrapping and magnetorotational instability (e.g., 49l-l5l| and references therein). It is 
worthwhile mentioning that such energetic MHD explosions are receiving great attention as a 



possible relevance to magnetars and collapsars (e.g., [52 



jwhich are presumably linked 
to the formation of long-duration gamma-ray bursts (e.g., 56| for a review). The grav- 
itational waveforms from the MHD explosions were found to have a quasi-monotonically 



57|, 



Such charac- 



increasing component that is produced by the bipolar outflows [48 , 
teristic waveforms with slower temporal evolution (< lOOHz) are likely to be detectable by 
future space-based detectors like DECIGO [59[, which are free from seismic noises at the 
lower frequency bands. 

Reliable modeling of rotating core collapse and the associated bounce signals should be 
done at least in three-dimension (3D), because nonaxisymmetric instabilities are expected 



to naturally develop in the postbounce supernova cores (e.g., [60|, l6l| for reviews). Probably 
the best understood type of instability is a classical dynamical bar-mode instability with a 
threshold of T/|iy| > 0.27 62l, l63|. Here T/|H^| represents the ratio of rotational kinetic 



energy to gravitational energy. In contrast to the high T/|14^| instability, recent work. 



some of which has been carried out in idea 
also in more self-consistent 3D simulations 



ized setups and assumptions |64|-|68J and later 



22 



23 



69| including microphysical EOS and a 



prescribed (Liedendorfer's) deleptonization scheme [70[, suggest that a differentially rotating 



We shortly call GW bursts emitted near bounce in the rotating iron cores as the bounce signals. 



PNS can become dynamically unstable at much lower T/|M^| as low as < 0.1. Despite 
clear numerical evidence for their existence, the physical origin of the low-T/|iy | instability 
remains unclear. It is most probable |66|, l7l| that the instabilities are associated with the 
existence of corotation points (where the pattern speed of the unstable modes matches 
the local angular velocity) inside the star and are thus likely to be a subclass of shear 
instabilitiesjj Even in the non-rotating case, it should be mentioned that t he p ostbounce 



core is no longer axisymmetric because of the growth of spiral SASI modes 
is considered to be a key ingredient for explaining the origin of pulsar's spin 



73 

73|. 



-|77|, which 



Since the vicinity of the PNS is an important cite for the non-axisymmetric instabilities, 
general relativity (GR) and accurate neutrino transport cannot be also neglected, not to 
mention 3D effects. In fact, Scheidegger et al.(2010) 23| found in their 3D post-Newtonian 
models that the GW amplitudes become 5-10 times bigger for models including deleptoniza- 
tion effects than those without. As is well known, the evolution and structure of the PNSs are 
very sensitive to the deleptonization episode via lepton/energy transport from inside (e.g.. 



78 



-|80|). In the above case, the deleptonization in the PNS results in the more compact and 



asymmetric core, leading to the stronger GW emission. A deeper gravitational well in GR 



simulations 



34 



81 



82| compared to the corresponding Newtonian models should affect the 



criteria of the low-T/|iy| instability. 

Despite the importance, recent multi-dimensional (multi-D) models both including GR 



and relevant microphysics are mostly 



models have been so far reported 



17 



37l 



imited to 2D (e.g., |2l|, yjsj), and only handful 3D 



8^. Ott et al. (2007) (l7| performed the first full 



3D GR simulations employing a realistic EOS and Liedendorfer's deleptonization scheme 



70| . In one of their 3D models with the largest GW amplitudes 17|, they reported the 



appearance of the low-T/|iy| instability after around tpb ~ 20 ms postbounce, which is 



characterized by the dominant azimuthal density mode oi m = 1 [6^, |65|- They found that 
the GW emission emitted in the direction of the spin axis become significantly stronger after 
the onset of the non-axisymmetric instability. One of the limitations in 17| is the use of 



the de^ 
phase 



eptonization scheme that was originally designed to be valid only in the prebounce 
70|. More recently, Ott et al. (2012) reported GR simulations employing a neutrino 



leakage scheme with 37| or without spatial symmetry assumptions in the computational 



Note that corotation resonance has been l ong known to the key ingredients in the accretion disk system, 
such as the Papaloizou-Pringle instabihty 



domain [8J] , nevertheless, detailed analysis of the GW signals (including models with rapid 
rotation) has not been reported yet. 

In this work, we study rotating core-collapse and the bounce GW signals by performing 
fully general relativistic 3D simulations with an approximate neutrino transport. The code 
is a marriage of an adaptive- mesh-refinement (AMR), conservative 3D GR MHD code devel- 



oped by [85|, and the approximate three-flavor neutrino transport code that was developed 



in our previous work (see 82| for details). The spacetime treatment in the code is based on 



86 



87|. Regarding 



the Baumgarte-Shapiro-Shibata-Nakamura (BSSN) formalism [see, e.g., 

neutrino transport, we solve the energy-independent set of radiation moments up to the first 

order and evaluate the second order momentum with an analytic variable Eddington factor 



(the so-called Ml closure scheme [88|). This part is based on the partial implementation of 
the Thome's momentum formalism, which was extended by [89'] in a more suitable manner 
applicable to the neutrino transport problem. By utilizing a nested grid infrastructure, an 
effective numerical resolution of our 3D models in the center is Ax ~ 450 m, which is as 



good as the most recent study by [37| and is better than our previous work [82] (Ax ~ 600 
m). By parametrically changing the initial angular momentum in the precollapse core of a 
15 Mq star ^], we follow the dynamics starting from the onset of gravitational collapse, 
through bounce, up to about 30-50 ms postbounce in this study. Albeit limited to the early 
postbounce phase (mainly due to the computational cost), we will show several interest- 
ing GW features, which come from non-axisymmetric spiral waves that develop under the 
influence of the low-T/jPF] instability and the spiral SASl. 

This paper is structured as follows. In Section HIl we briefly summarize the numeri- 
cal schemes, initial models, and the method how to extract the gravitational waveforms. 
The main results are presented in Section lllll We summarize our results and discuss their 
implications in Section HVl 

II. NUMERICAL METHOD 

This section starts with a brief summary about the basic equations and the numerical 
schemes of GR radiation hydrodynamics. We then move on to describe the initial condi- 
tions in Section m Bl and the method to extract the gravitational waveforms in Section fll Ct 
respectively. 



A. GR Radiation Hydrodynamics 



1. Metric equations 



The GR radiation hydro dynamic code developed in our previous work [82] consists of 
the following three parts, in which the evolution equations of metric, hydrodynamics, and 
neutrino radiation are solved, respectively (see, 82| for more details). Each of them is solved 
in an operator-splitting manner, but the system evolves self-consistently as a whole satisfying 
the Hamiltonian and momentum constraints. Note in this section that geometrized units are 
used (i.e. both the speed of light and the gravitational constant are set to unity: G = c= 1). 
Greek indices run from to 3, Latin indices from 1 to 3. 



Regarding the metric evolution, the spatial metric 7jj (in the standard (3+1) form: ds^ = 
—a'^dt'^ +'yij{dx'^ + f3''dt){dx^ + (3^ dt) , with a and /3* being the lapse and shift, respectively) and 
its extrinsic curvature Kij are evolved using the BSSN formulation 86|,|87|. The fundamental 
variables are 



— ln[det(7ij)] 



7ij 



A,, 



e-^n^. - ^1^JK) 



I 1. 



(1) 

(2) 
(3) 
(4) 
(5) 



The Einstein equation gives rise to the evolution equations for the BSSN variables as, 



{dt - Cp)^ij = 


-2aAij 


(dt - Cp)<P = 


-\.K 


{dt - Cp)A, = 


g-40 


a{Ri, - 87r7i^7,vT('^^t,i) 




-DiDjaf^ + a{KAij - 2 An 



(6) 
(7) 



zM 



(8) 



{dt - Cp)K = -Aa + a{AijA'' + K^S) 



(dt - /3''dk)r = levrf ^7,^ri,Tf 



(9) 



(total) 

+ ^r/3^. + /3^r^-2i^^a,-, (10) 

where £/? is the Lie derivative with respect to f3\ T'/^otai) ^^ ^^^ total stress-energy tensor (fluid 
+ radiation which we shall discuss in the next subsection), D denotes covariant derivative 
operator associated with •jij, A = D*Di,"trf" denotes the trace-free operator jJt,^ = (— a,0) 
is the time-like unit vector normal to the t = constant time slices. Following [91i], the gauge 
is specifled by the 1+log lapse, 

dta = I3'd,a - 2aK, (11) 

and by the Gamma-driver-shift, 

dt(3' = kdtt\ (12) 

here we chose k = 1. 



2. Radiation Hydrodynamics 

The total stress-energy tensor T^^^^i^ appeared in Equations (8-10) is expressed as 

^ (total) ~ ^ (fluid) "t" Z^ ^ (u) y^"^) 



where T^^^. and Tf f is the stress-energy tensor of fluid and neutrino radiation field, respec- 
tively. Note in the above equation, summation is taken for all species of neutrinos (z/g, z/g, i^x) 
with u^ representing heavy-lepton neutrinos (i.e. t'^, z/^ and their anti-particles). 
Starting from the definition of T^J), 



T^.f^ ^ E^^^n-n^ + F^^-^n^ + Ff^^n" + Pj^, 



(14) 



the evolution equations o: 
momentum formalism 92 1 



radiation energy {E(y^) and radiation flux (P")) in the truncated 
can be expressed as |89|, 






(15) 



and 






(16) 



where Q'^ denotes the source terms. In order to simplify the neutrino-matter interactions, a 
methodology of neutrino leakage scheme is partly employed at this moment (see [82| for the 
explicit expressions). 



By adopting the Ml closure J88|, the radiation pressure is expressed as 



p, io - '^ \pii 



(^) 



3(1 - X) „«j 



thin 



-* thick' 



(17) 



where x represents the variable Eddington factor, P^i^ and -Pthj^.^ correspond to the radiation 
pressure in the optically thin and thick limit, respectively. They are written in terms of Ey 



and F,, 



. For the variable Eddington factor x, we employ the one proposed by [88| 



X 



3 + 4F2 



5 + 2 V4 - 3F2 ■ 
E^ ' 



(18) 
(19) 



The hydro dynamic equations are written in a conservative form as, 

dtp, + diipy) = 0, (20) 

dtSi + dj{SiV^ + ae^'^PSi) = 

-Sodia + Skd.p'' + 2ae''^Sldi(P 
-ae^^{S,k - P7,fe)'9.7''V2 - e^^aQ^'^,^, 

(21) 
dtf + d,{Sov' + e^^P(t;^ + (3^) - p,v') = (22) 

-S^D'a + e'^^aQ^n^, (23) 

dt{p,Y,) + d,{p,Y,v') = pj:,, (24) 

where X = e^'^X, p, = pWe^'^, Si = phWui and Sq = phW"^ — P. p is the rest mass density, 
u^ is the 4-velocity of fluid, h = 1 +e + P/p is the specific enthalpy, f * = u" /v!'^ f = Sq — p*. 
Ye is the electron fraction, e and P is the internal energy and pressure, respectively. We 
employ the Shen EOS [93] for baryonic matter which is based on the relativistic mean field 
theory (see [82] for more details on the EOS implementation). 

B. Initial Models 



We employ a widely used progenitor of a ISMq star ((90|, model "sl5s7b2") and impose 
precollapse rotation manually to study its effect during the collapse, bounce, and the early 
postbounce phases in a controlled fashion. We assume a cylindrical rotation profile, 

u%^ = zuo^{no-^), (25) 

where u^ = ^u\ + m^ is rotational velocity and Wq set here to be lO^cm is reconciled with 
results from stellar evolution calculations suggesting uniform rotation in the precollapse core. 
The initial angular velocity at the origin Qq is treated as a free parameter and we compute 
four models by varying fio = 0, 7r/6, '7r/2 and vr (rad s~^). We enumerate models as Rn, 
where n ranges from (non-rotating) to 3 (rapidly rotating), that corresponds to the four 
choices of the precollapse central velocity (see Table [T|). Bearing in mind that it is very hard 
to go beyond ID computations in stellar evolution calculations and their outcomes may 

10 



not be the final answer, the central angular velocity for a 15Mq star is predicted to range 



from 0.15 |9^ to 3 (rad/s) [95 



depending on whether the (prescribed) angular momentum 
transport is included or not [96|. The central angular velocity for models Rl and R3 is 
adjusted to be closely in the same range with the progenitor models (e.g., Table H]). We 
add a random 1% initial density perturbation for all the models. By doing so, we hope to 
see convective activities shortly after bounce in our 3D-GR simulations including neutrino 
transport (albeit approximate) that are currently able to follow at most ~50 ms postbounce 
for each model (limited by the currently available supercomputing power at our hand). 

The outer boundary is at the radius of 7500 km and nested boxes with 9 refinement 
level are embedded in the computational domain without any spatial symmetry. Each box 
contains 128^ cells and the minimum grid size at the origin is thus Ax = 450m which is the 
same resolution with [37| and better than our previous study with Ax = 600m ^8^]. In the 
vicinity of standing shock front R ~ 120(240) km, our resolution achieves Ax = 1.8(3.6) km 
which is a factor of 2 (4) coarser than 371] . 



Gravitational Wave Extraction 



From GR simulations, GWs are extractable either by gauge- dependent or -independent 
methods. 63| compared the waveforms emitted from oscillating neutron stars by the two 
methods and reported the quadrupole formula produces the waveform in a sufficient accuracy 
compared to the gauge-invariant method. They showed the quadrupole formula underesti- 
mates the wave amplitude by ~ 10%, but the phase is not changed drastically. More detailed 
comparison is recently reported in |97|, which also supports that the quadrupole approxima- 
tion is adequate for stellar collapse spacetimes with a PNS formation. Accordingly in this 



23 



63, 



98 



99( 1 for extracting GWs 



work we choose the conventional quadrupole formula 
from our simulations. 

In the qu adrupole formula, the transverse and the trace-free gravitational field /ij, is 



expressed by 



lOrupoie 
IO0I.I1OI 



D 



(26) 



In Eq. (!26|) . A+/x represent amplitude of orthogonally polarized wave components, e+/x 
denote unit polarization tensors and D is the source distance. Following |22|, we adopt the 



11 



same expressions for the wave amplitude A+/x as, 



A, 



21, 



TT 



(27) 
(28) 



where the quadrupole moment in the sp heric al coordinates leg, Icj,^, and Ig^j) are expressed 



into the Cartesian components lij as 22|, Il02 | 



Ile^ = (III cos^ + F^ sin^ + 2/JJ sin cos 0) cos^ 9 



YTT 



TT . 



+rj sin' 6 - 2{I'J cos + /i' sin 0) sin ^ cos ^ 



Jjj = Jjj sin' + 47 cos' : 

rTT _ / YTT YTT 



yz 



iZ^ sin 20, 



xy 



(29) 
(30) 



(4^ -4x ) cos ^ sin COS. 



+/JJ' cos ^(cos' — sin' i 



+iJJ' sin 6* sin — ly^ sin 6* cos i 



(31) 



with superscripts TT denoting projection into the transverse-traceless gauge. The quadrupole 
moment I^j and its time derivative are given as 63| . 



and 



1.. = ^ 



G 



/_/ Jf JU iXj \-Ai JU « 



kj = — I p*{v'x^ + x'v^)d X. 



(32) 



(33) 



The second time derivative lij is directly computed by taking the time derivative of Eq. (133 
numerically. 



According to |25|, l35|], we estimate GWs from anisotropic neutrino radiation A'^^,^ as 
A!;(e, t) = ^f dt' I dn'il + s{9')ci^')s{0 + c(^')c(0) 



/\12 



X 



[s{e')ci<i>')ciO - cie')siO]' - [sie')si<i>')] 



and 



A'AU) 



'X^h^iu)' 



2G 



(34) 



C ./o 



X- 



dt' J dQ'ii + s{e')c{<p')s{0 + c{e')c{0) 

s{9')s{cP')[s{e')ci^')c{0 - ci9')siO] 
A9')c{<P')c{0 - c{e')s{0? + [s{9')sm^ 



x^l.-^h'- 



(35) 
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Here s{A) = sin(A), c{A) = cos(A) and FT-.' is radial energy flux of each neutrino species 
(z/ = h'e,i^e,^x) estimated at the radius of {Rl^ , 9' , (j)') . ^ denotes observer angle relative to 
the rotational axis and ^ = and 7i/2 for polar and equatorial observers, respectively. 

III. RESULTS 

A. Hydrodynamic features 

To summarize hydrodynamic features of our models, we first show evolution of the central 
(rest-mass) density and rotational /3 parameter (i.e., the ratio of rotational kinetic energy 
(T) to gravitational potential energy (W) near bounce in Fig{T] (see also table |T] in which 



several key quantities are summarized). Note that T and W is respectively defined as [103 1 



T=^ f p.hu^v'^dx^, (36) 

where v^ = u^/u^ and 

W = Mbar - MadM + ^int + ^rad, (37) 

with Mbar, Mabm, -Ejnt, and -Erad being the total baryon mass, ADM mass, internal energy, 
and neutrino radiation energy, respectiveljo. The top panel of Fig.l shows that the central 
densities for all the models at bounce exceed nuclear density (~ 2.8 x 10^*^ gcm~^). As 
consistent with p/7|, 137|, multiple bounce does not appear also in this study, which was 
otherwise often observed in previous Newtonian simulations with simplified setups. From 
the bottom panel of Fig.l, the /3 parameter at bounce reach 0.15, 1.3 and 4.9 % for models 
Rl, R2, and R3, which is consistent with J37| (hereafter Ott+12). Ott+12 reported full 3D 
GR simulations with relatively similar schemes as ours, i.e. the BSSN formalism including 
neutrino leakage (but without neutrino heating). Since our model series of RO - R3 employ 
similar initial central angular velocity as their "sl2WH07j(0-3)" in order, it is suitable to 
make comparisons. For instance, their /3 parameters for "sl2WH07j(l-3)" at bounce are 
0.4, 1.6, and 5.1 %, which are in good agree men t with our counterpart models (R1-R3). 



^ Following an analogy of electro-magnetic energy 
radiation field, i?iad, as 



103j, we evaluated the contribution from neutrino 



{E^^^W ^ Fl)U,)^dx\ (38) 
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I 
-o2 



u 



n=7T (R3) 

n=7T/2 (R2) 

n=7T/6 (Ri) 
n=o (Ro) 




10 , ^ 20 
Tpb (n^s) 



FIG. 1. Time evolution of the maximum rest mass density (top) and the rotational (3 parameter 
(bottom) as a function of postbounce time (Tpb) for models RO(black line), Rl(green line), R2(red 
line), and R3(blue line), respectively. 

On the other hand, there exist some differences in numerical setups between our study and 
Ott+12, including the adopted EOS, progenitor models, initial perturbations, and spatial 
restrictions. Several important impacts on the GW emission will be mentioned later. 

Fig. Elshows several snapshots of 3D distribution of specific entropy s(A;b baryon^^) inside 
the central cube of 700km^ for model R3, the most rapidly rotating model in this study. 



TABLE I. Model summary. The second column represents the initial central angular velocity (J^ini) 
followed by Pmax.b (the third column) and /3b (the fourth column), which represents the maximum 
rest mass density and rotational /3 parameter at bounce, respectively. 



Model 
RO 
Rl 
R2 
R3 



r^ini (rad s ^] 


it/6 
tt/2 

TT 



b(10i^g cm~3) 
3.54 
3.52 
3.41 
3.28 



/3b 
2.3 X 10"^ 
1.5 X 10-3 
1.3 X 10-2 
4.9 X 10-2 
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FIG. 2. Several snapshots of entropy distributions {k^ baryon~^) in the central cube of 700km^ 
for model R3 (top left; tpb = 9.9 ms, top right; tpb = 19.9 ms, bottom left; ipb = 29.9 ms and 
bottom right; tpb = 48.9 ms). The contours on the cross sections in the x = (back right), y = 
(back left), and z = (bottom) planes are, respectively projected on the sidewalls of the graphs 
to visualize 3D structures. 



After bounce, the shock rapidly expands in the direction toward the equatorial plane and it 
then stalls at an (angle-averaged) radius of ~260km at Tpb ~57ms. For the rest of models, 
the shapes of the shock surface remain nearly spherical (for model RO) or mildly oblate 
(for models Rl and R2) and the average shock radius roughly stays ~150 km within the 

15 




FIG. 3. Same as Fig. [2] but for the final time snapshots for models RO (top left), Rl (top right), 
and R2 (bottom). 



simulation time Tpb < 30ms (Fig. [3]). These features can be also seen in Fig. H] which shows 
evolution of average (dash-dotted), maximum (thick), and minimum (thin) shock radii for 
all the models. Pushed by strong centrifugal forces, the maximum shock extent becomes 
largest for the most rapidly spinning model (model R3), which is consistent with Ott+12. 

To see the effects of rotation on neutrino emission, we present in Fig. [5] the neutrino 
luminosities and the average neutrino energies £^,i- Since we do not transfer the number 
density of neutrinos in the present scheme (see |82| for more details), e^^^i can be evaluated 



16 




Tpb (rins) 

FIG. 4. Evolution of angle-average (dash-dotted), maximum (thick) and minimum (thin) shock 
radii as a function of postbounce time (Tpb). Black, green, red, and aqua lines are for models RO, 
Rl, R2, and R3, respectively. 

only by the following approximate way. We first project the positions of the neutrino sphere 



defined in the Cartesian grids to the spherical polar grids as R^ 



for each neutrino 



z/g, t'e; i^x- Then we estimate Ss^^i by the matter temperature at the neutrino 



species of i 

sphere assuming that neutrinos stream freely outwards with possessing the information of 

the last scattering surface. Accordingly Sg^ in the spherical coordinates (i?, 6', 0) is expressed 



as 



^s^,i{R, 0, 4>) = £u,i{Ru, 



i\" 1 Y )i " ^ 



Here e^^j in the right-hand-side denotes the neutrino energy at Ru{ 



(39) 



which is estimated 



by 



keT 



F3{Vu,0) 



(40) 



where Fk{ri^,f3) is the k-th Fermi-Dirac integral, r]^ = ^y/ksT is the degeneracy parameter 
with /ij, and T representing the neutrino chemical potential and matter temperature, re- 
spectively on the neutrino sphere, and /3 is the relativistic factor. The isotropic equivalent 
neutrino luminosity (e.g., 28|) is then evaluated by AttR'^FJ'. toward arbitrary polar angles 
for a given radial direction. 

As is expected, the neutrino luminosity and the average neutrino energy show little 
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FIG. 5. Postbounce evolution of the average neutrino energies Es^^i (top) and the isotropic equiv- 
alent neutrino luminosities (bottom) for models RO (left) and R3 (right). Thin and thick lines are 
for an observer along polar (z axis) and equatorial (x axis) directions. Black, red, and aqua lines 
represent electron, anti-electron, and heavy-lepton neutrinos, respectively. 

observer-angle variations for the non-rotating model RO and the peak z/g luminosity {L,yJ 
reaches ~ 3.6 x 10^^ erg s~^, which is quite similar to the non-rotating model in Ott+12. The 
neutrino energies of each neutrino flavor yield to the standard hierarchy (i.e. e^^ < Sp^ < e^^) 
within our simulation time. As seen from the right panels of Figure 5, the neutrino luminosity 
and the average neutrino energy for model R3, on the other hand, show a clear directional 
dependence. The peak Ve luminosity (L^,^ ~ 4 x 10^'^ erg s~^) toward the polar direction is 
approximately 10 % higher, compared to that along the equatorial direction (bottom right 
panel). By comparing with the luminosity for the non- rotating model RO (bottom left panel), 
the luminosity for model R3 generally becomes higher (lower) toward the polar (equatorial) 
direction. Due to the competition, the total luminosity becomes slightly smaller weaker 
(~ 6%) for model R3 compared to the non-rotating counterpart. By comparing the average 
neutrino energies (top left and right panels in Fig.5), the difference between each neutrino 
species becomes more remarkable for model R3 (i.e., the average neutrino energy becomes 
higher (lower) along polar (equatorial) direction, while preserving the mentioned energy 
hierarchy). All of these features are predominantly because of the rotational flattening of the 
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FIG. 6. (Left) Gravitational waveforms from matter motions for equatorial and polar directions 
(left-half and right-half, respectively) with two polarizations (top-half (x) and bottom half (+)) 
for all the computed models. Color represents the initial rotation rates; blue (model R3), red (R2), 
green (Rl), and black (RO), respectively. Note that the vertical scale is different for the bottom 
left panel. (Right) Same as the left panels but for the gravitational waveforms from anisotropic 
neutrino emission (sum of the contributions from all the neutrino species). Note that, to show 
a clear contrast between matter and neutrino signals, the wave amplitudes from neutrinos are 
multiplied by a factor of 10. 

core, by which the neutrino spheres of all the neutrino flavors are formed more deeper inward 
(outward) along the polar axis (equato r), preferentially enhancing the neutrino emission 
along the polar direction (e.g., 104j - |l06| ). In addition, the Uf, lightcurve near at the epoch 



of neutronization becomes more broader for model R3, which reflects the longer dynamical 
timescales (tdyn) for niodels with larger initial angular momentum|j. These features are in 
good agreement with those obtained in Ott+12[37|. 



B. Gravitational- Wave Signatures 

Now we are ready to discuss GW signatures in this section. After we shortly summarize 
the overall waveform trend, we perform detailed analysis on several new GW features that 
come genuinely from non-axisymmetric motions from the subsequent sections. 



^ This is because the central density (pc) supported by the centrifugal forces becomes smaller (idyn oc pc ) 
for rapidly spinning models (e.g., top panel of Fig. 1). 
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The gravitational waveforms from matter motions for all the computed models are shown 
in the left-half panels of Figure |6l Note in the panel that 

A+/xI = A+/x(0 = O, = 0) (41) 



and 



A+/xII = A+/X {9 = n/2, = 0) (42) 

(e.g., Eqs. (12 71128 p and f l34H35|) ) represent the quadrupole wave amplitudes with two polar- 
izations (x, +) for equatorial (denoted as "I" in the following) and polar (as "11") directions, 
respectively. 

From the lower left panel, typical GW features of the so-called type I waveforms (e.g., 
[l9|) are clearly seen in the A^ll mode, that is, a first positive peak just before bounce 
precedes the deep negative signal at bounce, which is followed by the subsequent ring-down 
phase. The wave amplitude for model R3 is in the range of —200 < A+II <100 cm near 
bounce. This is again comparable with that in the counterpart model "sl2WH07j3" in 
Ott+12. From Ott+07 who employed the same Shen EOS (but with a more simplified 
treatment for deleptonization) , their counterpart model with the /3 parameter ~ 6.7% at 
bounce has —240 < A^ll <90 cm, which is also in good agreement with ours. 

The right-half panels in Figj6] show the waveforms from anisotropic neutrin o em ission 



(Eqns. (IM|531) ). Due to the memory effect inherent to the neutrino GWs (e.g., [l07l |). the 
waveforms exhibit much slower temporal variation compared to the matter GWs. From 
the lower left panel, the largest amplitude of the neutrino GWs is obtained for the A+II 
mode in model R3 (blue line), which reaches ~ 6 cm with a quasi-monotonically increasing 
trend during the simulation time. This characteristic feature was already reported so far 
either in 2D simulations employing detailed neutrino transport |27| or in 3D simulations 



with idealized setups 35|. This comes from the (mentioned) stronger neutrino emission 
along the rotation axis due to the deformed neutrino sphere in the rapidly spinning cores. 
Except for model R3, the wave amplitudes from neutrinos stay around a few cm during our 
short simulation time. 
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C. Non-axisymmetric imprints 

Now let us have a look at again the matter GW signals for the A_^./^l modes (emitted 
toward the pole, e.g., the right-half panel of Figure 6 {left)). Note that these signals come 
genuinely from non-axisymmetric matter motions. The blue lines (model R3) in these panels 
show an oscillating behavior with two different modulation timescales (iTmod), firstly in a 
very short timescale (rmod ^ 1 ms) promptly after bounce within tpb < 8ms, and secondly 
in a relatively longer modulation {r^od ^ 10 ms) after that. These features are only clearly 
visible for our rapidly rotating model R3 (faintly visible for our moderately rotating model 
R2 (red line)). 

For model R3, the wave amplitudes in the first epoch (tpb < 8ms) reach |y4+/xl| ~20 
cm, while non- and slowly rotating models (RO and Rl) meanwhile produce very little GW 
emission {\A1\ < 1 cm). By comparing with the luminosity curves in Fig.5, this epoch is 
shown to closely correlate with the neutronization phase. During this epoch, the prompt 
shock propagates rapidly outward with capturing electrons and dissociating infalling iron 
group nuclei until the prompt shock stalls at around tpb ~ 10 ms. As seen from Fig. 4, 
the shape of the shock for all the models is rather spherical at this time, regardless of 
the difference in the initial rotation rates. In such a short duration after bounce, possible 
reason of producing the non-axisymmetry is less likely to be the low-T/|14^| instability, not 
to mention the SASI. We consider that the nonaxisymmetric matter mot ions come from 



the random perturbations imposed in the initial density profile (e.g., [2J, Il08j ). Because 
of the angular momentum conservation in the collapsing prebounce core, non-axisymmetric 
components seeded in the initial random perturbations have larger angular momentum with 
time. This leads to the sizable GW emission shortly after bounce before the GWs from 
prompt convection becomes more dominant. Fig. [7] shows the (absolute) GW amplitude for 
model RO with (black line) or without the initial seed perturbations (red line). As can be 
seen, the model without the perturbations (red line) shows much smaller growth rate of the 
wave amplitudes {\A\ ~ 0.01 cm) near around the neutronization phase (tpb < 8 ms) before 
prompt conve ction becomes active later on. Due to paucity of multi-D stellar evolution 
models (e.g., |l09| ). we still do not really know whether an 1% level of the precollapse 
density perturbation imposed in this study is realistic or not. More systematic study is 
clearly needed to unambiguously identify how much information we can extract about the 
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FIG. 7. The absolute GW amplitude of \A1\ (toward the polar direction, Eq. ()42p ) for non-rotating 
model RO with (black line) and without (red line) initial seed perturbations. 
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apse inhomogeneities (such as about their amplitudes, spatial fluctuation lengths, and 
) from the non-axisymmetric GWs. We also think that a correlation analysis with the 
neutrinos signals at neutronization should be complementarily helpful in this context. 

As the neutronization phase comes to an end at around tpb ~ 8 ms (Fig. 5), the prompt 
shock simultaneously stagnates, triggering the entropy-driven prompt convection behind. 
After that, the waveforms for moderately and rapidly spinning models (R2 and R3) exhibit 
a longer modulation with slightly larger amplitudes {\A\ ~ 20 cm) compared to those near 
the neutronization phase. Non- (RO) and slowly (Rl) rotating models emit roughly ~ 10 
times stronger GWs. These amplitudes during prompt convection are in accord with previous 
studies, \A\ ~ 15 cm in [StI, Is^, \A\ ~ 20 cm in [3^] and \A\ ~ 8 cm in 23]. 

As seen from blue line in Fig. 6 (left panel), the waveform of model R3 has a clear sinusoidal 
modulation, which possesses a tt/2 phase shift between A^l and Ay.1. This feature persists 
until the end of simulations (tpb ^ 50 ms in this study). To understand the origin of 
the sinusoidal signature, we present a spectrogram analysis in Fig. [81 In the figure, the 
angle-dependent characteristic strain /ichar 



h 



char I, 



F) 



2 G I dE{e,, 



(43) 



^ Unfortunately (exactly by the same reason for multi-D progenitor models), this is beyond the scope of 
the present work to run a wide variety of initial models due to their computational cost. 
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is plotted (e.g., [l6|, |32|, |36|), where D represents the source distance that we assume as 
D = 10 kpc (unless otherwise stated), F is the GW frequency, and dE{9, (f))/dF is the GW 
spectral energy density. dE{6, (j))/dF is given b}|I 

dE irc^ ^ 



dF 4G 



A+iF^ + \A,iF)\' 



Here A{F) denotes the Fourier component of A{t) 32| that is defined as 



A+/.{n 



A+,^{t)H{t-T)e 



-2mFt 



dt, 



(44) 



(45) 



where the width of the Hann window r is set as 10 ms. 

Equator 




10 20 

Tpb (tis) 



10 20 

f pb ("is) 



FIG. 8. Color-coded GW spectrograms (logio(^char)i Eq. (j43]) ) for models R3(top) and RO{bottom). 
Left and right panels are for polar [9, cj)) = (0, 0) and equatorial {9, cp) = (vr/2, 0) directions, 
respectively. The model name with the polarized mode is given in the lower left corner of each 
panel (for example, "R3 (AI)" in the top left panel). The contour drawn by black solid line in each 



11C| 



panel corresponds to a threshold beyond which the signal-to-noise ratio (SNR) for KAGRA 
exceeds unity (for a Galactic event). To guide the eyes (see text for details), the horizontal dotted 
line (200 Hz) is plotted. To clearly present the postbounce GW signatures, we set the maximum 
value of the color scale as logio(/ichar) = —21, while the wave amplitude is actually more higher 
near bounce (logio(/'-char) ~ —20.3 in the All mode (top right panel)). 



^ The angle dependence (6*, </>) is omitted for the sake of brevity. 
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FIG. 9. Same as Fig. [8]but for models R2 {top) and Rl (bottom). 

The spectrogram for the equatorial GW (top right panel) near bounce (tpb ~ ms) 
shows a power excess (colored by yellow to reddish regions) in the range between 



:00 < 

F < 1000 Hz, which is associated with the best-studied type I signals (e.g., [17|, ll8|, |20|)c. 

Not surprisingly, this feature is hardly seen in the non-rotating model (bottom right panel 

in Fig.8). Since the deviation from spherical symmetry for model RO is very small in the 

early postbounce phase, little angular variations are seen in the GW spectrogram (compare 

bottom left with bottom right panel). On the other hand, the spectrograms for model R3 

have much clearer angular dependencies (compare top left with top right panel). The polar 

GW spectrogram (top left panel) has two distinct power-excess islands that are enclosed 

by black solid line (satisfying the signal-to-noise ratio: SNR> 1 for a Galactic source by 

the second generation detectors) . The typical G W frequency of the first island is about ~ 

1000 Hz between < tpb < 10 ms (a reddish zone near the top left corner in the top left 

panel), and it is around 200 — 250 Hz for the second island with its excess clearly visible 

from tpb > 10 ms (e.g., a horizontal reddish stripe covering its peak around 200 — 250 Hz in 

the panel). 

The first island, which is narrower (in width and height) than the second one in the 

frequency-time domain, originates from the initial seed perturbations that we mentioned 

^ In line with previous works, these type I signals are likely within the detection limits of the next generation 

detectors for a Galactic source, when the /3 parameter /? at bounce exceeds ~ 1% (e.g., model R2 in this 

work) . 
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above. It is interesting to point out that the resulting GW amphtudes satisfy SNR> 1 
(enclosed by black solid line) for a Galactic source. In addition to the well-studied equatorial 
GWqj (associated with the type I signals), our results suggest that the polar GWs just 
after bounce have also a unique signature, which is produced by the percent levels of the 
precollapse density fluctuations seeded in the rapidly rotating cores. 

Regarding the second island, the characteristic frequency Fchar = 200 ~ 250 Hz for model 
R3 (tpb > 10 ms, top left panel)!^ is higher than those (Fchar ~ 100 Hz) in models RO and 
Rl (compare Figs. 8 and [9]). Note that the lower characteristic frequency (Fchar ~ 100 Hz) 
observed in our non-rotating and only mildly rotating models is consistent with recent results 



by [3J| who performed 2D GR simulations including detailed neutrino transport. According 
to their analysis, the GWs during prompt convection are predominantly generated by radially 
propagating acoustic waves above the PNS. With the typical sound velocity there (Cg ~ 10^ 
cm s~^) and the shock radius (-Rshock ~ 100 km), the typical frequency in our 3D-full-GR 
results (-Fchar ~ C'sZ-Rshock ~ 100 Hz) can be also reasonably estimated. On the other hand, 
the characteristic frequency for the second island (-Fchar = 200 ~ 250 Hz, most clearly visible 
for model R3, e.g., top left panel in Fig. 8) is shown to be systematically higher than the GW 
component solely from the propagating acoustic waves. From the next section, we look into 
the reason of the higher frequency in more detail. Touching on the detectability, the relevant 
GW frequencies (regardless of the difference) are in the range of -Fchar ~ 100 — 250 Hz, which 
are close to the maximum detector sensitivity for the second generation interferometers (e.g.. 



KAGRA and Advanced LIGO [6|,l8|). 

D. GW Emission from One-Armed Spiral Waves 

In this section, we are going to discuss that the characteristic frequency (-Fchar = 200 ~ 
250 Hz); clearly visible in the GW spectrogram for model R3) has a tight correlation with 
the one-armed spiral waves that develop under the influence of the SASI and the low-T/|W^| 
instability. We will also show that these features can be naturally explained by the acoustic 
feedback between the stalled shock and the rot ating PN S surfaces, which is thus reconciled 



with the well-established picture of the SASI 



m|-[ll3|) 



^ We shortly call GWs seen from the polar (equatorial) direction as polar (equatorial) GWs, respectively. 
^° e.g., the horizontal dotted line that closely divides the island into two. 
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FIG. 10. Similar to the right panels of FiglHl but contributions to the first time derivative of 
the mass quadrupole moment lij are presented for different radial locations from the innermost 
{R < 20km, top), 20km< R < 60km (middle), to the outer region (R > 60km, bottom), respectively. 
Note in this panel that GWs seen from polar direction are plotted. 

Similar to the right panels in FigJHl but we plot in Fig. [10] contributions to the first time 
derivative of the mass quadrupole moment (Eq.(l33l)) from different radial locations (seen 
from the polar direction). The top, middle, and bottom panels represent /ichar evaluated 
by the spatial integral in the following range, within i? < 20 km, 20 < i? < 60 km, and 
i? > 60 km, respectively. It can be seen that the GW emission in the second island (again, 
the horizontal (Fchar = 200 ~ 250 Hz) reddish zone after around Tpb ~ 10 ms) are radiated 
mainly from above the PNS [R > 60 km, bottom panel) during prompt convection. Note 
that the spatial location itself is similar to the one reported in [3j] for their non-rotating 
model. 

To see how the source of the strong GW emitter evolves with time, Fig{TT] shows the 
integrand ip in the quadrupole formula (top panel); 



V' = (A+I)2 + (AxI)' = (/..-/, 



yyj 



{21 



xy) 



(46) 



and the specific entropy (bottom panel) at selected time slices. From the upper panels 
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FIG. 11. Colormaps of the integrand tp in the quadrupole formula (top) and entropy (per baryon, 
bottom) on the equatorial plane at selected time slices (denoted in the top right corner in each mini 
panel) for model R3. In the upper panels, the central region i? < 20 km is excised to show a clear 
contrast. Note that the color scale is in a logarithmic scale and it is normalized by the maximum 
value of 4>. 

(regarding ip), it can be seen that the one-armed spiral wave starts to form at around (e.g., 
Tpb = 6.74 ms, top middle panel (upper)) for model R3, and then it keeps rotating in a 
counter-clockwise manner (compare Tpb = 26.75 ms and Tpb = 28.75 ms) until the end of 
simulations. From the two panels (at Tpb = 26.75 and 28.75 ms), the spiral wave rotates 
around ~ 90° during 2 ms, thus the rotational period (frequency) can be estimated as 
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FIG. 12. (iyg/it) Space-time diagrams of the one-armed spiral components of vorticity (top) and 
acoustic waves (bottom) (see text for more details). {Right) To guide the eyes, dotted arrows are 
inserted to illustrate the ad vective- acoustic cycle on top of the contour lines from the left panels. 

^rot ~ 8 ms and 2/T^ot ~ 250 Hzlll|. Note that this is close to the mentioned GW frequency 
-^char = 200 ~ 250 Hz (e.g., Fig. [8]). Comparing with the lower panels, the spiral mode is 
seen to connect between the standing shock front (the reddish high-entropy regions between 
the outer blueish regions and inside) and the vicinity of the PNS (~ 50 — 60 km in radius). 
In the next section, we are moving on to seek the possible origins of the spiral waves in more 
detail. 

E. Possible Origins of One-armed Spiral Waves 

By carefully looking at the top panels (Fig. 11), kinks are formed at the triple points 
where the standing shock and the spiral wave meets. Such type of morphology has been 
ubiquit ously obs erved in previous simulations aiming to unravel the nature of the SASI 



73l . uSi, Ill3l-lll8l| ). To check more carefully in our results whether the advective-acoustic 
cycle 119|, [l20| is running or not, we plot in Fig{T2]the space-time diagram of vorticity (top 



Here 2 in the numerator comes from two polarized wave modes (+/x' 
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panel, denoted by "Advective" ) and acoustic amplitudes (bottom panel, by "Acoustic" ) for 
model R3. Here the two quantities are evaluated as 

aiw) = -^ , (47) 



M^ 



where X = -ojV ■ (f 0e^)o |l20| and X = p {p the pressure) are for the vorticity and acoustic 



amplitudes, respectively, ro = ya?+y^, = tan~^ {y / x) , and Yi^rn{0-i4>) is the spherical 
harmonics. Note that purely sloshing SASI modes were studied in detail based on 2D 
simulations 1201 ]. but here we primarily focus on the (/,'m) = (1, 1) mode to extract the 



characteristic pattern of the one-armed spiral wave. 

To guide the eyes, the right panel in Fig|T2] illustrates propagation of the advective (red 
lines) and acoustic waves (blue lines), respectively. Just behind the shock (see, point "A" in 
the top right panel), the generated vorticities go down to the PNS surface (e.g., close to the 
point "B"). This is supported by the arranged direction (e.g., from top-left to down-right 
direction in the panel) of the bumpy islands, each of which is separated by black contour 
lines. The reddish regions at i? ~ 10 — 20km after Tpb ~ 20 ms (see the bottom left 
colormap) represent strong generation of the acoustic waves near the PNS. They propagate 
outward until they hit the shock. This is seen from the direction of the greenish stripes in 
the colormap (bottom left), which is symbolically drawn by blue lines in the bottom right 
panel. As indicated by red lines after point A (top right), the adve ctive- acou stic cycle is in 



operation subsequently. Above features are in good agreement with 1191 . Il20| (see references 



Jierein), although the S ASI signa tures are not as clearly discernible as in the previous 2D 



12 1| or 3D simulations 1131 . IllSl ] due to the short postbounce evolution in this study. 

In order to more clearly specify the non-axisymmetic structures, we monitor density 
profiles at a given radius w in the equatorial plane, which are decomposed into its azimuthal 
Fourier components as 

J T p{'^A,z = 0)e^ 
jl'^ p{w,(t),z = 0)d(l) 

For a mildly rotating model (Rl, left panel of Figure 13), the Cartesian m = 4 background 
noise (blue lines) dominates over the other modes almost always inside the PNS [zo = 20 
km, thick line) and at the shock {w = 130 km, thin line). On top of the Cartesian noise, the 



Am = ^^^£^-^^^ ^:7^-^. (48) 



Here e^ is a unit azimuthal vector 
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FIG. 13. Time evolution of the normalized amplitudes of the density for different azimuthal modes 
{m = 1,2,3, and 4) for models Rl (left) and R3 (right). Thick and thin lines are extracted at 
■CD = 20 and 130 km, respectively. 

linear growth of m = 1 mode (inside the PNS for model R3, thick blue line) is clearly seen in 
the right panel, which gradually transits to the saturation phase as the prompt convection 
phase sets in (tpb > 10 ms). Note that from the thin blue line (right panel, for w = 130 km), 
the presence of the one-armed spiral waves behind the shock (compare Fig. [TT]) is evident. 
The dominant m = 1 mode inside the PNS is considered to be of a remarkable signature of 



the low-/3 instability (e.g., 17|, |23[). Our results show that not only the (mentioned) spiral 



SASI but also the low-/3 instability coexist in the rapidly rotating supernova cores. 

Back to FigJT2l stronger acoustic waves are generated at Tpb ~27, 35, and 44 ms near the 
PNS approximately in the time interval of ~8 ms. This timescale is close to the characteristic 
GW frequency of the spiral waves (-Fchar = 200 ~ 250 Hz. e.g., in Fig. 6). In order to better 
address the origin of Fchar ~ 200-250 Hz, FigjMl shows angle-averaged frequencies associated 
with the rotational velocity (i^rot, dotted lines) and the acoustic wave (r2rot + ^aco, solid 
lines) in the equatorial plane. Note that they are respectively defined as 



n 



rot 



2vr/ 






y 



Q =2- 

2t[^ x^ + y^ 



(49) 
(50) 



and the sum (firot + ^aco) is a measure to estimate the propagation timescale of the acoustic 
wave on top of the rotating medium. As seen, the acoustic frequency (solid lines) shows 
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FIG. 14. Profiles of fiaco + ^rot (solid lines) and r^rot (dashed lines) along the equatorial direction 
(x axis) at different time slices (see text for more details). 
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^200-250 Hz in the range of 60 km < /? < 120 km, which is actually in good 



agreement with the narrow-band GW emission seen in Figs. (TOlandlTTl In the regions above 
the PNS (60 km ^ R ^ 120 km), it it also shown that firot ~ 100 Hz is as high as fiaco- The 
frequency difference (~ 100 Hz) between the Doppler-shifted acoustic frequency (firot + ^aco) 
and the purely acoustic one (fiaco) can naturally explain the phase shift regarding the peak 
GW frequency between model R3 and the remaining models with smaller initial angular 
momentum. 

Above results also show that the purely acoustic frequency Jlaco ~ 100 Hz is less sensitive 
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to the initial angular momentum. In fact, the GW spectra in the literature [23 
generally peak at around ~ 100 Hz during prompt convection. We thus speculate that 
significantly higher GW peako, if observed in the spectrogram during prompt convection, 
might be a possible signature of rapid rotation, which is unobservable if not for the GW 
astronomy. 



^^ To make it possible, a correlation analysis with neutrino signals (e.g., 37L l4ll . ll22l |l should be indispensable 
to specify the epoch of core-bounce. 
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F. Detectability 



To discuss detectability, FigJTSl shows the charac terist ic GW spectra /ichar of our selected 
models with the design noise curves of initial LIGO 122| . Advanced LIGO |6|, and KAGRA 



110| ). assuming a source distance of 10 kpc. 
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FIG. 15. Characteristic GW amplitudes /ichar of our selected models ( R3 (in green), R2 (in red). 



and RO (in black)) compared with the strain sensitivity of initial LIGO 122], Advanced LIGO a], 
and KAGRAJllOt] at a source location of 10 kpc. Left and right panels are for a spectator along 
polar and equatorial directions, respectively. Near the bottom left corner, the initial angular 
velocity of each model is provided for a reference. Note that due to the small contribution from 
the neutrino-originated GWs, only the matter contribution is plotted here. 

For our rapidly rotating models (e.g., R2 and R3), /ichar along the equatorial direction 
(right panel) is generally within the detection limits of the advanced detectors for an assumed 
distance of 10 kpc. The corresponding signal-to-noise ratio (SNR) is approximately greater 
than ~ 10 over a wide frequency range 100 ^ F < 1000 Hz. In accord with previous 2D 
and 3D results 






171 . Il9l . IST] , faster initial rotation (but not too much) increase the chance for 



detection. The spectral peak appears at Fpcak ~ 620 Hz for models R3 and R2 (green and 
red lines in the right panel), which is associated with the maximum spike seen in the type I 
waveforms near bounce (e.g., Fig|6] and the spectrogram in Fig 18]). Quantitatively this is in 
good agreement with Ott+07 [l7| who obtained -Fpeak ~600-700 Hz for their counterpart 3D 
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model, in which the same Shen EOS was employed (but with much idealized microphysical 
treatment). 

Seen from the polar direction, the GW spectra have two distinct peaks (e.g., green and 
red lines in the left panel of Figure fT5|) . The first peak, appearing at -Fpcak ~ 700 Hz (model 
R2) and ~ 1000 Hz (model R3), is emitted during the neutronization phase (tpb < 10 
ms). As mentioned earlier, this component totally vanishes in axisymmetry and it reflects 
the presence of precollapse density inhomogeneities. The second peaks are seen around 
-^peak ~ 100 Hz for uou- to moderately-rotatiug models (red and black lines in the left 
panel). They are predominantly determined by the characteristic timescale of the acoustic 
waves traveling between the PNS and the stalled shock during prompt convection (see section 
unci) . Due to rapid rotation, the acoustic frequency is shifted upward Fpcak ~ 200 Hz in 
model R3 (green line in the left panel). As mentioned in the previous subsection, this is due 
to the Doppler effect (e.g., Eqns. (I49l[50l) ) of the acoustic waves emitted on top of the spiral 
waves. The SNR of the first peak (-Fpeak = IkHz) for a Galactic source is at most SNR~ 10 
even for model R3, whose sensitivity is limited by shot noise at high frequencies. On the 
other hand, the SNR of the second peak achieves as high as SNR~ 100 (for model R3), 
because it is close to the maximum detector sensitivity for the advanced interferometers. 

IV. SUMMARY AND DISCUSSIONS 

We have studied properties of GWs from rotating core-collapse of a ISMq star by per- 
forming 3D-full-GR hydrodynamic simulations with an approximate neutrino transport. By 
parametrically changing the precollapse angular momentum, we paid particular attention to 
the effects of rotation on the GW signatures in the early postbounce evolution. Regarding 
neutrino transport, we solved the energy-independent set of radiation energy and momen- 
tum based on the Thome's momentum formalism. In addition to the matter GW signals, 
we took into account GWs from anisotropic neutrino emission. In addition to common 
GW signatures obtained in previous 2D axisymmetric studies, our results showed several 
non-axisymmetric features in the waveforms which can be explored only by 3D simulations. 

Among the common features are the type I waveforms emitted along the equatorial 
direction, which have comparable amplitudes with those in previous 3D GR studies [l7|, |37| ■ 
The wave amplitude reaches a few xlOO cm for our most rapidly rotating model R3. For 
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a Galactic source, this is well within the detection limits (with the SNR > 10) over a wide 
frequency range (100 ^ F < 1000 Hz) of KAGRA or Advanced LIGO. in our moderately to 
rapidly rotating models. The peak GW frequency (Fpeak ~ 620 Hz) in the GW spectra is 
also comparable to Ott+07 who employed the same Shen EOS for hadronic matter. During 
prompt convection, the gravitational waveforms do not show any qualitative differences 
except for the most rapidly rotating model R3. The wave amplitude is \A\ ~ 10 cm (non- 
and slowly rotating models) and 1^41 ~ 20 cm (moderately one) which are consistent with 
previous GR simulations |3J, |37| and these signals reach the SNR~ 8, independently of the 
observer direction. 

We also studied neutrino luminosities, the average neutrino energies, and the waveforms 
associated with anisotropic neutrino emission. Rotation makes the neutrino luminosity 
and the average energies higher toward the polar direction, and conversely lower along the 
equatorial direction. These features are consequences of rotational flattening of the central 
core and neutrino spheres. The stronger neutrino emission toward the rotational axis leads 
;o a quasi-monotonically increasing trend in the wave amplitude, as qualitatively similar to 



271] who employed more sophisticated neutrino transport. The impact of rapid rotation on 
the neutrino GWs was found only for the polarized wave mode A^ll in our most rapidly 
rotating model R3. 

Our findings clearly show that non-axisymmetric instabilities play an essential role in 
determining the GW signatures in the rotating postbounce evolution. By analyzing the GW 
spectrograms, the GW emission toward the rotational axis has two distinct features, which 
appears at the higher {F ~ 700 — 1000 Hz) and lower frequency domain {F ~ 100 — 200 
Hz), respectively. The higher one near bounce, which appears clearly in our moderately 
and rapidly rotating models, comes from the precollapse density inhomogeneities. For a 
Galactic source, their signals achieve SNR~ 5 — 10 for the advanced detectors. The lower 
one seen in non- to moderately rotating models is originated from spiral waves that develop 
under the control of the spiral SASI and the low-T/|iy| instability. In our most rapidly 
rotating model, the lower peak frequency is shifted upward to F ~ 200 Hz which is due to 
the Doppler effect (e.g., Eqns. ( H9l l50l) ) of the acoustic waves emitted on top of the spiral 
waves. Regarding the detectability, the relevant GW frequencies are close to the maximum 
detector sensitivity for the second generation interferometers (e.g., KAGRA and Advanced 
LIGO), which is thus expected to be detectable (SNR~ 100) for a Galactic source. 
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Finally we would like to discuss some of the limitations of the present work. For a 
more quantitative GW prediction, the simplifications in the employed neutrino reaction 
(and the gray transport scheme) should be improved, which we regard as the most urgent 
task (Kuroda, Takiwaki, and Kotake in preparation). We need to conduct a convergence 
check in which a numerical gridding is changed in a parametric manner, although it is too 
computational expensive to do so for our 3D-GR models at present. In addition, we need 
to run as many models as possible to study the dependencies of the precollapse density 
inhomogeneities, the progenitor dependence, the initial rotation rates and magnetic fields 
on the GW signals. Some good news is that we have access to the "K-computer" which is 
among the world- fastest Peta-scale supercomputers. By utilizing it, we hope to study these 
important themes one by one in the near future. 
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